A unifying definition of synchronization for dynamical systems 
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We propose a unifying definition for synchronization. By 
example, we show that the synchronization phenomena dis- 
cussed in the dynamical systems literature fits within the 
framework of this definition. 



Synchronization between dynamical systems has been 
an active research topic since the time of Huygens. It is 
a phenomenon of interest to fields ranging from celestial 
mechanics to laser physics, and from communication to 
neuroscience 0. 

Over the last decade, a number of new types of syn- 
chronization have appeared: chaotic synchronization 
phase synchronization ||, lag synchronization 0, and 
generalized synchronization Q, to mention only a few. 
This is in addition to the classic examples of synchroniza- 
tion in periodic systems [^|J^]. Many of these have been 
experimentally observed in a single system Q . Synchro- 
nization is often categorized on the basis of whether the 
coupling mechanism is uni-directional or bi-directional. 
Stable synchronization with uni-directional coupling has 
been called synchronization by an external force (for fre- 
quency synchronization) and master-slave synchroniza- 
tion (Pecora and Carroll Q). (It has recently been shown 
that, if the synchronized systems are identical then there 
is no essential difference between uni-directional and bi- 
directional synchronization |9|.) 

Although there have been several attempts jl(],[nj , no 
successful definition of synchronization currently exists. 
The definition in use is an ever increasing enumerated 
list. When a "new type" of synchronization arises, its 
name is added to the list. We believe that "definition by 
example" is an untidy situation which should be replaced 
by a single definition that encompasses all of the known 
examples. 

In this letter we propose a unified definition which ac- 
counts for all types of synchronization between finite di- 
mensional systems. Although we explicitly discusses syn- 
chronization between two continuous time dynamical sys- 
tems, our results can be extended to N continuous time, 
or N discrete time systems. Therefore, our results apply 
to a larger class of phenomena than the one we explicitly 
discussed. 

To construct this definition, assume that a large sta- 
tionary deterministic dynamical system is divided into 



two sub-systems, 



dx 



(i) 



Here, x G IR dl and y S lR d2 are vectors that may have 
different dimensions. The phase space and vector field of 
the large system is formed (in a natural way) from the 
product of the two smaller phase spaces and vector fields. 
Examples of phenomena that can be described by Eq. (|l|) 
are ubiquitous. 

Colloquially, synchronization means correlated in-time 
behavior between different processes. In fact, the Ox- 
ford Advanced dictionary |[12[ , defines synchronization 
as "to agree in time" and "to happen at the same time" . 
From this intuitive definition we propose that synchro- 
nization requires: (1) Separating the dynamics of a large 
dynamical system into the dynamics of sub-systems. (2) 
A method for measuring properties of the sub-systems. 
(3) A method for comparing the properties of the sub- 
systems. (4) A criteria for determining whether the prop- 
erties agree in time. If they agree then the systems are 
synchronized. The remainder of this letter formalizes 
this intuitive definition of synchronization by explicitly 
addressing each requirement, and applying the proposed 
definition to examples. 

We begin by separating the dynamics of the large sys- 
tems into the dynamics of sub-systems. Let 4>{zq) de- 
note a trajectory of the large dynamical system, given by 
Eq. ([]]), with initial condition, z = [x ,y ] G M 1 ®lR d2 . 
Respectively, curves <P x (zq) and <p y (zo) are obtained 
from this trajectory by projecting away the y and x com- 
ponents. We say that <p x (zo) and (j> y {z ) are "trajecto- 
ries" of the first and second sub-systems of Eq. (|l|). In 
this context we have separated the trajectories <p x (zq) 
and <p y (zo) from </>(zo)j rather than constructing 4>(zq) 
from 4> x {zq) and <}> y {z ). 

To discuss measuring properties of the sub-systems let 
X denote the space of all trajectories of the first sub- 
system, and consider a mapping g x : X ® 1R — > M k . 
The first M represents time, and is included so that g x 
can make explicit reference to time. We say that the 
mapping, g x , is a property of the first sub-system. The 
image of [cf) x (z ),t] G X®IR under the mapping g x is the 
result of measuring the property of the first sub-system. 
It will be denoted by g(x) G JR k . Similar definitions 



can be made for the second sub-system. The following 
examples make these notions less abstract. 

For synchronization, a property of a sub-system that 
is often of interest is frequency. Measuring the property 
means calculating a numerical value for the frequency. 
Hence, u> x = g(x). Other properties of interest are the 
coordinates of a sub-system at time t. Measuring the 
properties means determining numerical values for the 
coordinates. Hence, x(t) — g(x). Experimentally, g x is 
the quantity being measured, and g(x) is the value of 
the measurement. These examples show that a property 
can be a long time average, or a quantity whose value 
depends implicitly on time. Furthermore, the dimension 
of the measurement, k, can take on different values de- 
pending on the property being measured. Notice that it is 
reasonable to say g x is a property of the first sub-system 
because g(x) is obtained without explicitly referring to 
any other sub-system. 

Finally, we discuss the notions of comparing the prop- 
erties, and determining when they agree in time. We say 
the function h : M k ® M k -> M k compares the measured 
properties of the two sub-systems, and the two measure- 
ments agree in time if and only if h[g{x), g(y)] = 0. Be- 
low, a norm is used to determine this last requirement. 

With these preliminaries in place, we offer the follow- 
ing definition for synchronization: 

Definition The sub-systems in Eq. (Qj are synchronized 
on the trajectory (j){zo), with respect to the properties, 
g x and g y , if there is a time independent function h : 
R k ®R k -> R k such that 

\\h[g(x),g(y)]\\=0, 

where || • || is some norm. 

A subsequent definition removes details of initial condi- 
tions and trjectories: The sub-systems are synchronized 
with respect to the properties g x and g y if the previous 
definition holds on all trajectories. The subsequent defi- 
nition is what many papers in the literature call synchro- 
nization However, as we and others have shown, syn- 
chronization depends strongly on the trajectory JUpjj. 
Therefore, the trajectory dependence in the first defini- 
tion can not be ignored. 

We claim that this definition naturally follows from 
the intuitive definition of synchronization, and that it 
encompasses all of the interesting examples found in the 
literature. A strength of this definition is that the proper- 
ties and comparison function are not specified, a priori. 
As shown below, different applications require different 
properties and comparison functions, and those that are 
suitable for one application are often completely unsuit- 
able for another. For example, the following comparison 
functions all appear in the literature 



h[g(x),g(y)] 



lim — 

T-oo T 



t+T 



[g(x(s))-g(y(s))]ds. (4) 



h\g(x),g(y)] = g(x) - g(y) 
h[g(x),g{y)] = lim [g(x) - g(y)\ 

t—>OG 



Similar breadth occurs with properties. The same two 
sub-systems may be synchronized with respect to some 
properties, yet not synchronized with respect to other 
properties. A lack of breadth is the Achilles heel of pre- 
vious definitions. Typically, they fail because they, a pri- 
ori, specify properties and/or comparison functions that 
must be applied to all types of synchronization. There- 
fore, it is easy to find examples where the selected prop- 
erty and/or comparison function is inappropriate. 

Some may be concerned that the definition is too gen- 
eral. We argue, via analogy, that generality is also 
strength. Perhaps the most useful concept in theoreti- 
cal physics is a vector space. The definition of a vector 
space is as general as the one proposed for synchroniza- 
tion Q] . The definition does not specify what constitutes 
a "vector" or the operation "+" . Thus, a range of things 
from matrices to Fourier series to bras and kets are vec- 
tors in their respective vector spaces. The definition only 
insists that the set of "vectors" obey a specific series of 
abstract rules. If a set obeys these rules then it is a vector 
space, and the considerable power one obtains from that 
knowledge can be employed. (Group theory is another 
example of an extremely useful concept in physics whose 
definition is abstract.) Our definition, gives an explicit 
list of four tasks and a condition that must be satisfied 
for synchronization. Like the definition of a vector space 
(or a group) it provides a structural framework that can 
be used for subsequent research. We submit that this 
structure is an improvement over the current situation of 
failed definitions and enumerated lists. 

The remainder of this letter demonstrates the utility 
of the definition by discussing well-known examples. 

Frequency Synchronization 

Sub-system properties used in frequency synchronization 
are frequencies. If the trajectory is mostly rotation about 
an axis then the measured frequencies (lo x = g(x) and 
Lu y = g(y)) are located at power spectra peaks associ- 
ated with the average rotation of the signal. Examples 
of such dynamics include, periodic motion, systems with 
phase coherent chaotic attractors (like Rosslcr 15 1), or 
systems with Silnikov dynamics |lrj[ ]. For these examples, 
phase modulation contributes weakly to the dynamics, 
and chaos (if it exists) results mainly from amplitude 
modulation. 

The measurement function is typically h[g(x),g(y)] = 
n x oj x — riyLUy (where n x and n y are integers). Synchro- 
nization implies that the frequencies of the sub-systems 
are commensurate 



n x to x - riyUy = 0. 



(5) 



(2) 
(3) 



Many text books discuss frequency synchronization for 
periodic systems M. 



Frequency synchronization between coupled chaotic 
systems with bi-stable attractors has also been exam- 
ined JlTj. (The Lorenz and double scroll attractors are 
examples of bi-stable attractors.) In Ref. [jl7), the prop- 
erties are the average frequency of switching between the 
two lobes of the attractors, and frequency synchroniza- 
tion on a trajectory occurs if Eq. ^ is satisfied. Our defi- 
nition works for this example. Many other definitions fail 
because, either the definition of phase is ambiguous jl8) , 
and/or \\x — y\\ need not remain small |10|. 

Another example is frequency synchronization between 
a chaotic bi-stable attractor and a periodi c sy stem. The 
communication method of Hayes et. al. |Lj| labels the 
lobes of the attractor as or I, and uses small control 
signals to produce a trajectory which encodes the mes- 
sage. The obvious choice for the sampling rate of the 
receiver is the mean switching frequency of the chaotic 
system. Therefore, a periodic receiver with frequency u> v 
is synchronized to a chaotic transmitter with switching 



frequency uj x if u> x 



This type of synchronization 



fits into our definition, but does not seem to fit into any 
previous definition. 

Frequency synchronization compares properties that 
are long time averages of the trajectory. Therefore, it 
is a loose restriction on the dynamics of the sub-systems. 
In particular, it does not restrict the instantaneous values 
of the coordinates x and y. All remaining examples com- 
pare properties whose measured values depend implicitly 
on time. 

Phase Synchronization 

Phase synchronization involves sub-system properties 
called "phases". If the dynamics is chaotic and phase 
coherent then one can introduce cylindrical coordinates, 
and unambiguously define the phase as the angle coordi- 
nate, 4>(t). However, other applications define the phase 
via a Hilbert transform, in which case the phase may not 
be uniquely defined on the sub-system Also, there 
are examples where the measured phase, <fi(t), is a vec- 
tor obtained from a trajectory using none of the previous 
methods @. 

If the measured properties are given by g(x) — <fi x (t) 
and g(y) = <p v (t) then the most common comparison 
function is 



h[g(x) )g (y)} = U[e,(g(x)-g(y))} 



(6) 



Here, U(u, v) is a vector with a-th component 
U a (u,v) — © [u a — \v a \], and is the unit step func- 
tion. Equation (|^) says that synchronization means 
\4>xa — 4>ya\ < e Q , so e is the maximum tolerable sep- 
aration between the components of the phase. The value 
of ||e|| is usually small, but can not be set a priori because 
its size is application dependent (2Cj . If "phase slips" oc- 
cur then a comparison function using a time average like 
that of Eq. (|j) is necessary || . 

Phase synchronization only compares the phase vari- 
ables. In the synchronous state the phases are locked, but 



the amplitudes can remain chaotic and relatively uncor- 
rected. Our definition includes phase synchronization. 
In contrast, definitions which focus on \\x — y\\ jl(J fail 
because phase synchronization does not restrict ampli- 
tudes. Likewise, any definition that forces one to use 
a specific type of phase will fail because phase is not 
uniquely defined. 

Identical Synchronization 

This is the most frequently discussed form of synchro- 
nization within the nonlinear dynamics community [Q. 
Here the sub-systems are identical, and the properties are 
the phase space variables, g(x) = x(t), and g(y) = yit). 
Most discussions in the literature use Eq. (^) as the com- 
parison function [p],^0[. 

However, if the dynamics of the sub-systems are 
chaotic then bursts (sudden loss and recovery of syn- 
chronous motion caused by unstable periodic orbits 
within the attractor) occur on chaotic trajectories for 
some forms of coupling p3]. Applications which can 
not tolerate bursts demand high quality synchronization, 
where Eqs. (^) is used as the comparison function [ BOj . 
If bursts are tolerable then a hybrid of Eqs. (jij) and (g) 
can be used. 

An engineering application called "dead-beat" syn- 
chronization (only possible in discrete time dynamical 
systems) uses Eq. (g) as the comparison function 
This type of synchronization is often used to describe 
systems whose measured properties are restricted to a 
finite symbolic alphabet p2j. 

In their seminal paper, Afraimovich, Verichev, and Ra- 
binovich [^3| generalized identical synchronization in two 
different ways. 

Lag Synchronization 

Two sub-systems are lag synchronized if their measured 
properties lag each other by a fixed amount of time, r. 
A trivial example is when the measured properties are 
g(x) — x(t) and g(y) = yit + r), and the compari- 
son function is Eq. (0). For this example, the second 
sub-system follows the same trajectory as the first sub- 
system, but is t units of time behind. 

A nontrivial example is Ref. Q. In this paper the 
measured properties are g(x) = Xi(t) and g(y) = yi(t + 
t) (the first components of x and y). The comparison 
function is h[g(x),g(y)] = K ([g(x) - g{y)} 2 ), where K 
is a constant and (•) is a time average. For this example, 
the sub-systems are not identical and S 2 (t) = \\h\\ = 
for a non-zero value of t. 

Lag synchronization also occurs if, instead of a con- 
stant value of r, one uses g(y) = y[T(t)], with T : M — > 
M a homeomorphism with lim^oo 2^1 = l pjjfl . 

The second generalization in Ref. p3j follows from the 
observation: If the sub-system are identical then the set 
x = y defines an invariant manifold in the phase space 
of the large system. 



Generalized Synchronization 

The literature is not consistent when discussing gener- 
alized synchronization. Most papers say that general- 
ized synchronization occurs if the measured properties 
are g(x) = x, g(y) = y, and the comparison function, 
h, is given by 



h[g(x),g(y)]=H[g(x)]-g(y) 



(7) 



where H is a smooth, invertible, time independent func- 
tion 11 1 . Roughly speaking, the sub-systems are gen- 
erally synchronized if y(t) = H[x(t)]. The equation, 
y = H(x), defines an invariant manifold in the phase 
space of the large system, and one can determine the 
state of one sub-systems from the state of the other sub- 
system [ p"H . 

However, Rulkov et. al. Jsj examined an example 
where the sub-systems have the same functional form 
but different parameter values. Their numerical and ex- 
perimental evidence indicates that it is possible to have 
stable frequency synchronization on a trajectory and not 
have generalized synchronization in the sense discussed 
above. For their example, one sub-systems oscillated 
twice for every oscillation of the other sub-system (i.e., 
oj\/u>2 = 2). This implies that it is impossible to con- 
struct a smooth invertible mapping y = H{x). There- 
fore, the definitions in Ref. fail. 

This example illustrates that the definition of general- 
ized synchronization needs to include J?'s with a finite 
(or perhaps countable) number of branches. Similar con- 
clusions arise from Ref. |24[ . Because we only require the 
functional reationship, h[g{x), g(y)] = 0, between the 
properties of the sub-system, both of these examples are 
naturally contained within our definition of synchroniza- 
tion. In particular, we do not insist that this relationship 
have the form g(y) — H[g(x)]. Therefore, we argue that 
the definition we have proposed is a more natural defini- 
tion for generalized synchronization. 

In conclusion, we have describe four tasks that are re- 
quired for synchronization. Based on this, we proposed 
a unified definition of synchronization between finite di- 
mensional dynamical systems. We claim that this defi- 
nition encompasses all examples of synchronization dis- 
cussed in the literature, and that it offers a common lan- 
guage and framework that can be used to discuss different 
types of synchronization. 
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